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Abstract 

Recently there have been experimental indications that solid 4 He might be a supersolid. We discuss the relation 
of supersolid behavior to ring exchange. The tunnelling frequencies for ring exchanges in quantum solids are 
calculated using Path Integral Monte Carlo by finding the free energy for making a path that begins with the 
atoms in one configuration and ends with a permutation of those positions. We find that the exchange frequencies 
m solid 4 He are described by a simple lattice model which does not show supersolid behavior. Thus, the PIMC 
calculations constrain the mechanism for the supersolid behavior. We also look at the characteristics of very long 
exchanges needed for macroscopic mass transport. 
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Recent torsional-oscillator observations by E. 
Kim and M. H. W. Chan on solid 4 He, both in the 
disordered absorbed vycor and in bulk He [1], 
have revived interest in the super-solid phase. A 
supersolid [2] is characterized by both long-range 
translational order and by a non-classical response 
to rotation (NCRI). In attempting to understand 
these results, one should first assume that the 
observed phenomena is described by equilibrium 
thermodynamics of pure bulk solid helium. For 
such a model, it is possible to perform numeri- 
cally exact Path Integral Monte Carlo simulations 
and explore its basic properties. If such properties 
agree with experiments, that will provide strong 
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support for the model. If one finds disagreement, 
then the basic assumption is invalid. The results of 
such a calculation have been recently published [3]. 
Here, we elaborate and update that publication. 

We consider the possibility that bulk solid hep 
helium, assumed to be free of defects, could have an 
equilibrium supersolid phase. Because of its light 
mass, solid helium has large quantum zero point 
fluctuations causing instantaneous ground state 
defects. At any instant of time, a good fraction 
of atoms are closer to a neighboring site than to 
their home site. However, the absence of an atom 
from a lattice site is not sufficient for having a su- 
persolid; if the empty site is always accompanied 
by doubly occupied site, there can be no mass cur- 
rent. However, Chester[4] proved that any Jastrow 
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(i.e. pair product) wavefunction of finite range 
has both BEC and vacancies. Although Jastrow 
wavefunctions do not make a good description of 
the ground state, one can multiply by a periodic 
trial function[5,6] and achieve a crystalline, Bose 
condensed state and a reasonable description of 
the energy, though not as low as a wavefunction 
without BEC. Though this wavefunction of the 
solid is not too bad for static properties, other 
properties, particularly off-diagonal long range 
order (ODLRO), equivalent to BEC, is not nec- 
essarily accurately given. There is no reason to 
think that even a perfect match of the diagonal 
properties constrains the off-diagonal properties. 

Path integrals give a rigorous way of calculat- 
ing on and off-diagonal properties. The partition 
function of N bosons is: 

Z = ]vTE / ' dR{R\e-f> H \PR) (1) 
• P J 

where H is the Hamiltonian, (3 the inverse tem- 
perature and R — {ri, r 2 , . . . rjy}. We will assume 
the helium atoms interact with a semi-empirical [7] 
pair potential. The helium potential is known to 
be accurate; many calculations agree with exper- 
imental measurements, without the adjustments 
needed in all other atomic and molecular systems. 
Using the PIMC method, superfluidity and freez- 
ing, happen naturally at the right density and tem- 
perature, without imposing them in any way as is 
done in a wavefunction based approach. 

To perform the PIMC calculations, the density 
matrix operator is expanded into a path, R(t) 
where t (0 < t < j3) is imaginary time. Then the 
superfluid density (i.e. the number of atoms not 
moving with the walls of the apparatus) is pro- 
portional to the mean squared winding number, 
p s oc< W 2 > where W = dt dvi(t)/dt 
is the winding number of the paths around the 
periodic boundary condition. Note that only ex- 
changes on the order of the sample size contribute 
to the superfluid density; local exchanges make 
no contribution. Prokof'ev and Svistunov[8] argue 
that long windings are impossible without vacan- 
cies and interstitials that are spatially separated. 

The technical complications of path integral 
calculations [9] concern ergodicity of the random 



walk, and finite size effects: one has to take the 
limit as N — > oo. If one does a PIMC simulation 
of a 48 atom (3x4x2) hep supercell in PBC, one 
finds a superfluid density of about 3% at melting 
density (molar volume 21.04 cm 3 ) and about 1.2% 
at 55 bars (molar volume 19.01 cm 3 ). However, 
simulations of a larger cell 180 atoms, show zero 
superfluid density. With current algorithms, it is 
difficult to be sure that the lack of winding paths 
has a physical origin, or is due to a lack of ergod- 
icity within the random walk. Winding number 
changes involve the permutation of a group of 
particles which stretch across the box. One needs 
to update simultaneously 5 atom coordinates to 
change the winding number in the 180 cell. The 
acceptance probability for such a collective move 
becomes very small. To avoid the problem of esti- 
mating rare events, we turn to a PIMC approach 
which directly estimates individual exchange 
probabilities. 

The Thouless[10] theory of exchange in quantum 
crystals assumes that at low temperatures, the sys- 
tem will almost always be near one of the iV! ar- 
rangements of particles to lattice sites. There will 
be rare, rapid, tunncllings from one arrangement 
to another in an amount of "imaginary time" 7 P . 
Hence, we can label the particles with their ini- 
tial lattice sites and drop the N\ in Eq. (1). The 
partition function is written as a sum over permu- 
tations of lattice sites onto themselves. We break 
up the permutation into cyclic exchanges P = 
{pi,P2 ■ ■ - Pn}- We can then write: 

np 

z = z Q J2l[f P M- (2) 

P i=l 

where Zq — J dR(R\e^ l3H \R) is an uninteresting 
phonon partition function at low temperature and 
the contribution for an elementary cycle is: 

W) = yJ dR(R\e- fJH \ P R) = J P (3. (3) 

The exchange is localized in imaginary time so 
the path integral is proportional to the number of 
crossing times: f p ((3) oc j3 ; the coefficient, J p is the 
exchange frequency. We have also assumed that in- 
dependent cycles do not interfere with each other, 
reasonable because the exchanges are rare and are 
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Table 1 

Calculated exchange frequencies in hep 4 He at a molar 
volume of 21.04 cc in units of fiK. The notation of the 
exchange is from Rogcr[13] as are the WKB actions for 
3 He. The column J 4 g are for /3 = 1, r = 0.0125 and N = 48 
and N = 180 for a larger cell. 7 p is the exchange time in 
units of K . 



p 


name 


^48 % 


error 


^180 % 


error 


7p 


WKB 


2 


nn 


1.7 


2 


3.74 


7 


0.32 


9.54 




nn' 


2.5 


6 


4.50 


5 


0.31 


8.60 


3 


T 


5.6 


4 


8.0 


5 


0.29 


7.33 






0.62 


5 


1.25 


6 


0.34 


8.83 




T' 


2.0 


4 


4.9 


5 


0.32 


7.75 


4 


Kp' 


0.74 


9 


0.59 


6 


0.19 


9.00 




Kp 


0.84 


8 


0.71 


6 


0.4 


9.38 




Ksq' 


1.25 


4 


2.2 


8 


0.34 


8.09 



unlikely to be close in both space and imaginary 
time. 

We have developed methods [11, 12] to calculate 
the exchange frequencies in quantum crystals, pri- 
marily for solid 3 He. In that system, ring exchanges 
of 2-6 atoms give rise to the magnetic properties 
at low temperatures. The PIMC calculations were 
important in giving a microscopic justification of 
the multiple exchange model. The PIMC values are 
in good agreement with experimental data on 3 Hc 
and are thus expected to be accurate in the closely 
related solid 4 He. 

In Table 1 are shown the exchange frequencies 
for 2, 3 and 4 atom exchanges in hep 4 He. For 
these calculations, we use an imaginary time step 
of 0.0125-fT -1 and a value of imaginary time of (3 = 
2K^ 1 (160 time slices). For most of the exchanges 
we find an exchange time j p = 0.3i^ _1 . We have 
studied systems with 48 and 180 atoms; the larger 
value is necessary for convergence of the exchange 
frequencies as is shown in the table. The values are 
quite small e. g. J 2 ~ 3[iK at melting density. We 
note that in an hep crystal, there are often sev- 
eral alternative ring exchanges of n atoms[13]. For 
example, one can have a 2 particle near-neighbor 
exchange in the basal plane n or out of the basal 
plane n'. Similarly a 3-body ring exchange can ei- 
ther enclose atoms in nearby planes T or not T*. 
The frequencies for these various alternatives can 
differ by an order of magnitude. 



The small cyclic exchanges in Table 1 are quite 
different from the long exchanges needed to get 
a supersolid. Accordingly, we have performed ex- 
change calculations of 50 different exchanges in- 
volving from 5 to 10 atoms. All exchanges are be- 
tween nearest neighbors, since calculations show 
that next-nearest exchanges are much less proba- 
ble. We obtain accuracies on the order of 5% for 
the 5-particle exchanges and 10% for the 9 par- 
ticle exchanges. More than half of the exchanges 
we studied involve winding around the cell bound- 
aries, important because they are representative of 
the type of exchanges in a supersolid. 

Fig. 1 shows the results of calculations of the fre- 
quency of the simplest winding exchanges: straight 
line exchanges in the basal plane. These calcula- 
tions were done in an hep supcrccll by varying the 
number of unit cells in the x-direction and keep- 
ing a fixed width in the y and z directions. We find 
that the exchange frequencies decrease exponen- 
tially with the length of the exchange with an ex- 
ponent of a = 2.64 near the melting volume 21.04 
cm 3 , and a — 3.14 at the volume 19.01 cm 3 corre- 
sponding to P <~ 60 bars. 

In order to understand whether the Bose- 
Einstein transition occurred in supcrfluid 4 He, 
Feynman introduced a lattice model [14]. The 
Feynman-Kikuchi (FK) partition function is: 

Z= Z ^e- a £r=i L « (4) 
p 

where the sum is over all closed non-intersecting 
polygons drawn on a lattice by connecting near- 
est neighbors. Calculations were done on a cubic 
3D lattice, first by approximately by Feynman and 
Kikuchi and then numerically by Elser[15]. It was 
found that the superfluid transition occurs as a be- 
comes smaller from a localized state (small cycles) 
to large cycles for a c « 1.44. To understand this 
critical value, consider the free energy of adding a 
link to an existing loop. Adding a link costs prob- 
ability e~ a but the entropy of the new link is 5 
in a cubic lattice; hence a c ~ ln(5) = 1.6. The 
exact critical value is smaller because of the 'self- 
avoiding' restriction within a cycle and on overlap- 
ping cycles. Using the same argument for the hep 
lattice gives a critical coupling of a c w 2.3. Our 
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Fig. 1. The exchange frequencies (J in K) versus exchange 
length L p for straight line exchanges in the basal plane that 
wind around the periodic cell. The inset shows the lattice 
sites in the basal plane for L p = 6; the red atoms show 
the atoms midway through the exchange. The circles with 
error bars are the ln(J) at two molar volumes; 21.04 cm 3 
(black, open circles, solid line) and 19.01 cm 3 (red, solid 
circles, dashed line). The lines are least squares fits. 

determined slope in figure 1 is larger than this but 
not by much! 

However, we need a more realistic model than 
the FK model to take into account more details of 
the geometry of the exchange than just the number 
of exchanging atoms; note for example the large 
differences between the 2 types of nearest neigh- 
bor 4-atom exchanges, a parallelogram K p and a 
square K sq in the Table. We assume that it is 
only the internal geometry of the exchange that 
matters [16] ; the detailed arrangement of the neigh- 
boring spectator atoms is much less important. In 
particular, we assume the the log of the exchange 
frequency is a sum over the internal angles 9k of 
the exchange: 

p 

J p = J exp[- ^2 a($k)]. (5) 
k=l 

In an hep lattice there are 7 possible angles be- 
tween two nearest neighbor displacement vectors 
(neglecting possible dependence on motion with 
respect to the basal plane), but 9 — only oc- 



Fig. 2. The fitted angular dependance of the action vs. 
cos 4 (0/2) at the melting density , 25 bars (open circles). 
The angle is the interior angle of a vertex as shown. The 
fits were done using 50 different exchanges in a 180 atom 
cell. The solid red circles correspond to a pressure 60 bars. 
The two lines show a critical supcrsolid Hamiltonian with 
(solid line) and without (dashed line) angular dependance. 

curs in the pair exchange which we neglect since it 
does not contribute to the superfluid density. We 
determine the values of the 7 parameters by fit- 
ting to the PIMC exchange frequencies J p . Good 
fits are obtained; the model predicts the exchanges 
frequencies with an accuracy of about 20%. (Note 
that the data base does not include the small ex- 
changes in table 1. The frequencies for small ex- 
changes are more sensitive to the neighbors.) There 
is a strong preference for exchanges that proceed 
in a straight line versus ones which have sharp an- 
gles; for them the incoming and outgoing particles 
are more likely to collide. We find that the fitting 
coefficients depend linearly on cos 4 (9/2) so that in 
the generalized F-K model we assume: 

Tip 

J p = J cxp[— aL p — a' cos i (9 i /2)\. (6) 
i=i 

With the improved model, we can return to the 
question of whether solid 4 He is a supersolid. Be- 
cause of the directionality of exchange, the prob- 
ability of retracing is small and we can ignore the 
self intersections, allowing us to find the probabil- 
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Fig. 3. Phase diagram of the generalized Feynman-Kikuchi 
model. The 2 points shown are for a density near melting 
and at high pressure. 

ity of long exchanges by formulating it as a Markov 
process: the diffusion of a particle in an hep lat- 
tice. Note that on the basal plane, exchanges can 
go straight (6* = 180), while an exchange out of 
the basal plane must zig-zag (8 = 146), reduc- 
ing the exchange probability, thus the diffusion is 
anisotropic. Then, an atomic displacement only de- 
pends on the previous atom's displacement, not on 
the position in the lattice. (Because the hep lattice 
has a basis, we have to label the displacements con- 
sistently on the A and B planes so that the tran- 
sition probabilities are independent of the plane.) 
Let the (un-normalizcd) probability of a given dis- 
placement vector after n steps be 11^ where i 
refers to one of the lattice directions. Then it will 
satisfy the master equation 

n j™ +1) =rnf ) e a(<lrJj> (7) 

i 

The probability for long exchanges will correspond 
to the maximum eigenvalue A of its stationary 
state: 

^2lI i e a V t -W = AIL;. (8) 

i 

By symmetry LL can only depend on whether the 
direction is in the basal plane, II i, or out of the 
basal plane, II2. Define the partial sum Dij = 
e~ a ( 8i ~ e i} where are either in the basal 
plane (1) or out of it (2), e.g. Dn is the sum of 
the probability of both vectors being in the basal 
plane). The eigenvalue A, is the solution to the sec- 
ular equation: 

GDu - A)(D 22 - A) - D\ 2 = 0. (9) 



Table 2 

Parameters for the generalized F-K exchange model, Eq. 
5 as determined by fitting to PIMC determined exchange 
frequencies. 



v (cm 3 ) 


P (bar) 


Jo(K ) 


a 


a> 


A 


21.04 


25 


7.2 


2.64 


4.8 


0.31 


19.01 


60 


8.0 


3.14 


4.8 


0.25 



\=(D 11 +D 22 )/2 + D 12 [l + y 2 } 1 / 2 (10) 

where y = (D 22 - D n )/(2D 12 ) w -0.023. Putting 
in the model parameters, we find A = 0.303±0.005. 
The probability of having an exchange of length n 
will equal a prefactor times A", hence, since A < 1, 
our PIMC results imply that solid He will have 
only localized exchange and thus cannot be a su- 
persolid. Since y is small, the diffusion in and out 
of the basal plane are similar. Self-intersections 
and the presence of other exchange cycles will fur- 
ther decrease the probability of long exchange cy- 
cles and increase the critical value of A. (To esti- 
mate the effect of self-intersections quantitatively, 
we performed random walks on the hep lattice and 
counted non-intersecting walks. This line is shown 
in Fig. 2.) As seen on Fig. 3 the estimated value of 
A is much less than what is needed to allow for a 
supersolid. In addition, there is no indication that 
long exchanges become prevalent near melting. 

One assumption made in the calculations is that 
long exchanges are correctly being treated by the 
method to calculate J. We find that the exchange 
time is roughly independent of the number of ex- 
changing atoms, implying that an exchange of L p 
atoms happens by all of the atoms simultaneously 
sliding over to their new positions. As L p increases, 
such a mechanism seems unlikely There would be 
much more phase space for the exchange path if 
different portions of a long path exchanged at dif- 
ferent times. 

To understand the mechanism for long ex- 
changes, we have done special simulations of 
them. Though it is not possible to calculate the 
frequency by the existing algorithm, it is straight- 
forward to sample long exchange paths. In partic- 
ular, we have studied the exchanges shown in fig. 
1, where L p atoms slide along the basal plane in 
the x-direction, where 10 < L p < 50. The bound- 
ary conditions in imaginary time are R(0) = Z 
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and R{j3) = PZ where P is a cyclic exchange of 
atoms along a straight line. (The cell has a fixed 
extent in the y and z-direction.) We do single 
particle updates on all the atoms using the third 
level bisection procedure. Because of the ergodic 
problem, the initial exchange will never decay, so 
we can study it at will. In particular, we are inter- 
ested in whether a long exchange maintains the 
instanton character referred to above, or whether 
it becomes delocalized in time. 

In order to quantify the exchange process, we 
define a reaction coordinate of each atom involved 
in the exchange as r\i (t) = — 1/2 + (xi(t) — Xi(0))/a 
where a is the nearest neighbor distance and Xi (t) 
is the x-coordinate of atom i at imaginary time t. 
The reaction coordinate is defined so that r\i (0) = 
— 1/2 and r)i((3) = 1/2. Fig. 4 is a gray scale image 
of the reaction coordinates for a typical path with 
L p = 40. One sees a fairly sharp division between 
the initial lattice sites (in black on the left) and the 
final lattice position (light color on the right). We 
define the exchange time ei for a given atom as the 
solution to rji(ei) — 0. In case of multiple solutions 
we find the solution where the neighbors atoms are 
also exchanging. The exchange times for the path 
in fig. 4 are shown as dots and connected with a 
solid line. 

The resulting curve of exchange times can be 
used to characterize the exchange process. Nearly 
vertical portions, represent instantons. In the fig- 
ure, 16 atoms (labelled on the y axis as 27-40 1-3 ) 
are such a segment and the 5 atoms (8-12) another. 
However, in long exchanges, such segments are bro- 
ken by an interstitial vacancy (i-v) formation, a 
subsequent pinning of the interstitial and diffusion 
of the vacancy. An example is atom 16 followed by 
the diffusion of the vacancy (atoms 17-26). The ex- 
change rarely proceeds by a single i-v event. The 
jagged appearance of the exchange curve is confir- 
mation that vacancies and interstitials are strongly 
bound and cannot be separated by more than a few 
lattice spacings. This is different than the creation 
of a single i-v pair that is discussed in Prokof'ev 
and B. Svistunov[8]. In the path represented in 
fig. 4, at least 3 i-v fluctuations and 4 or more in- 
stantons are present. By forcing the long exchange, 
we create multiple i-v pairs as in a recently pro- 
posed supersolid mechanism [18]. However, such i- 
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Fig. 4. Gray scale image of the reaction coordinate of a 
winding path containing 40 atoms, as a function of particle 
number (y axis) and imaginary time (x-axis). The darker 
region on the left side have coordinates closer to the initial 
lattice sites, the lighter region on the right side is closer 
to the permuted lattice. The black dots show the reaction 
time: when each atom is midway between its initial and 
final positions. A vertical portion of the reaction time curve 
(such as at the top) represents a group of atoms that move 
as a unit. Horizontal (or jagged) parts of the curve represent 
interstitials. 




Fig. 5. The "roughness" of the exchange time curve ver- 
sus the number of exchanging atoms. See the text for the 
definition of roughness. The solid curve is a calculation of 
a simple Brownian bridge model of the roughness. 

v pairs arc not low energy excitations; otherwise 
they would have been detected in specific heat or 
other thermodynamic measurements[19]. 

To quantify this process, we examine the "rough- 
ness" of the exchange curve versus the system size. 
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Fig. 6. The auto-correlation (see text) of the exchange time 
vs separation of atoms. Shown are the results for L p =10 
(black), 20(red), 30(blue) and 40(green). The curve goes 
negative for large cycles because it must average to zero 
and has periodicity L p . The dots are analytic results for a 
Brownian bridge. 

We characterize the roughness by first sorting the 
exchange times and discarding the largest and 
smallest L p /10 values to make the statistics more 
robust. The roughness is then t r = (e max ~ e m i n ). 
The roughness grows as L]J 2 , see fig. 5. Also 
shown is the same measure of roughness com- 
puted on a Brownian bridge with L p elements 
(random ID diffusion returning to the start af- 
ter L p steps). Figure 6 shows the autocorrelation 
function c(k — j) = {SekSej) /t% of the exchange 
time for systems with 10 < L p < 40. The auto- 
correlation function is very close to that given by 
that of a Brownian bridge: c(x) = 6x 2 — 6x + 1. 
This is another indication that the exchange time 
is a random diffusion process. Hence, for long ex- 

1/2 

changes, the exchange time should grow as L p , 
because of the creation of a i-v pairs every 5 to 10 
atoms. The decay rate of the exchange frequency 
will still be exponential: J = exp(— aL p ) since 
such a law results from both an instanton process 
and the creation of multiple i-v pairs. 

In summary, PIMC-computed exchange fre- 
quencies for hep solid 4 He show that only localized 
exchanges will be present and thus should not ex- 
hibit the property of nonclassical rotational iner- 
tia. We find no indication that long cycles behave 
differently than short cycles. In a separate study, 
we modelled 4 He absorbed in Vycor[17] and found 
that a persistent liquid layer picture could explain 
the experimental findings. To explain the bulk 
4 He experiment one must look for other mecha- 
nisms, either non-equilibrium dynamical effects or 



more complicated lattice defects not allowed by 
the boundary conditions in our simulations. 
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